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Strong intergalactic shocks are a natural consequence of structure formation in the universe. These shocks 
are expected to deposit large fractions of their thermal energy in relativistic electrons (£ e — 0.05 according 
to supernova remnant observations) and magnetic fields (£b — 0.01 according to cluster halo observations). 
We discuss the expected synchrotron emission from such shocks, and the observational consequences for next 
generation radio telescopes such as the Square Kilometer Array. We present an analytical model, calibrated 
and verified based on a hydrodynamical ACDM simulation. The resulting signal composes a large fraction 
(up to a few 10%) of the extragalactic radio background below 500 MHz. The associated angular fluctuations, 
e.g. STi > 260({ e f B /5 x 10 _4 )(^/100 MHz)~ 3 K for multipoles 400 < I < 2000, dominate the radio sky for 
frequencies v < 10 GHz and angular scales 1' < 9 < 1° (after a modest removal of point sources), provided that 
£e£s > 3 x 10 -4 . The fluctuating signal is most pronounced for v < 500 MHz, dominating the sky there even 
for £ e £s = 5 x 10~ 5 . We find the signal easily observable by next generation telescopes such as the SKA, and 
marginally observable with present telescopes. The signal may be identified using cross-correlations with tracers 
of large-scale structure (e.g. 7-ray emission from intergalactic shocks). Detection of the signal will provide the 
first identification of intergalactic shocks and of the warm-hot intergalactic medium (believed to contain most of 
the baryons in the low redshift universe), and gauge the unknown intergalactic magnetic field. 



1. Introduction 

The gravitational formation of structure in the 
universe inevitably produced strong, collisionless 
shocks in the intergalactic medium (IGM), owing 
to the convergence of large-scale flows. In these 
shocks, electrons are believed to be Fermi acceler- 
ated up to highly relativistic (> 10 TeV) energies, 
limited by inverse-Compton cooling off cosmic mi- 
crowave background (CMB) photons [23] ■ The re- 
sulting 7-ray emission traces the large-scale struc- 
ture of the universe. Rich galaxy clusters should 
be detected by future 7-ray missions as 7-ray 
sources 13 66 70 in the form of accretion rings 
with bright spots at the intersections with galaxy 
filaments [221 - Recently, a possible association 
of 7-rays (measured by the Energetic Gamma- 
Ray Experiment Telescope, EGRET) with Abell 
cluster locations was identified at a 3cr confidence 
level (OJ, but see also [fTT]). 
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In addition to the inverse-Compton emission 
from intergalactic shocks, synchrotron radiation 
should also be emitted by the relativistic elec- 
trons, as they gyrate in the shock-induced mag- 
netic fields. The resulting radio signature is ex- 
pected to trace the structure of the universe at 
low redshifts (z < 1), and to be dominated by 
rich, young galaxy clusters. Indeed, extended ra- 
dio emission with no optical counterpart is ob- 
served in about 10% of the rich galaxy clusters 
(radio halos and radio relics, see |24|). and in 
more than a third of the young, massive clusters 
(with X-ray luminosity Lx > 10 45 erg s _1 , see 
|22|). The radio emission has been identified as 
synchrotron radiation from relativistic electrons, 
but there are different models for the origin of the 
electrons and the magnetic fields involved (for a 
recent discussion, see 0]). Observations of nine 
radio relics were used (20] to suggest an associa- 
tion between these sources and structure forma- 
tion shocks, focusing on the possibility that the 
relics are revived fossil radio cocoons originating 
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from nearby radio galaxies, re-energized through 
diffusive shock acceleration. Recently, large-scale 
diffuse radio emission was discovered around a 
filament of galaxies, possibly tracing an accretion 
shock on this scale jj|]. 

Waxman & Loeb (2000, [70]) have proposed 
a simple model, which allows one to estimate 
both the radio and the 7-ray signatures of in- 
tergalactic shocks, produced by electrons accel- 
erated from the inflowing plasma. Their model 
allows one to estimate the radio and the 7-ray 
backgrounds and their anisotropy characteristics, 
as well as the signature of individual clusters. The 
model uses dimensional analysis arguments to es- 
timate the properties of the virialization accre- 
tion shock of a halo, as a function of redshift z 
and halo mass M. Halo abundance estimates at 
different redshifts [such as the Press & Schechter 
(1974, halo mass function] , may then be used 
to calculate various observables. The model ap- 
proximates the strong accretion shocks as being 
spherically symmetric, and neglects weak merger 
shocks. A fraction £ e ~ 0.05 of the shock ther- 
mal energy is assumed to be carried by relativis- 
tic electrons, based on observations of supernovae 
remnants (SNRs, for a discussion, see (35] ). Ob- 
servations of > 0.1 /xG magnetic fields in cluster 
halos (for a recent review, see jUj) require that a 
fraction £g — 0.01 of the shock thermal energy be 
transferred into downstream magnetic fields. The 
model has been generalized for a ACDM universe 
by Keshet, Waxman & Loeb (2004a, EH]). The 
free parameters of the model were calibrated ac- 
cording to global features of a ACDM simulation, 
and verified to give radio and 7-ray signatures 
consistent with the signatures found directly from 
the simulation. Shock asymmetry was incorpo- 
rated into the model by introducing a geomet- 
rical correction factor, approximately measured 
from the simulation. 

Identification of radio or 7-ray emission from 
intergalactic shocks holds a great promise for ad- 
vancing current knowledge on shock formation 
in the IGM. ft should provide the first direct 
evidence for such shocks, revealing the underly- 
ing large-scale cosmological flows. When com- 
bined with 7-ray detection, the radio signal will 
provide a direct measure of the unknown mag- 



netic fields in the IGM, possibly shedding light 
on the processes leading to IGM magnetization. 
Emission from large-scale shocks traces the unde- 
tected warm-hot IGM at temperatures 10 5 K < 
T < 10 7 K that is believed to contain most of 
the baryons in the low redshift universe (see e.g. 
|17|). Moreover, the signal may be used to study 
non-thermal physical processes in the intergalac- 
tic environment, such as Fermi acceleration in low 
density shocks. 

In this chapter we discuss the radio signature of 
intergalactic shocks, and the observational conse- 
quences for next generation radio telescopes such 
as the Square Kilometer Array. In we review 
the model for radio emission from intergalactic 
shocks, and test its validity. In S|2]we discuss the 
observational consequences of the predicted radio 
signal, with emphasis for the SKA, and exam- 
ine confusion with various foreground and back- 
ground signals. Finally, g3l summarizes our re- 
sults and addresses their potential implications. 
We discuss methods by which the signal may be 
identified, and present some consequences of a fu- 
ture positive detection of the signal. 



2. Model 

Here we study the extragalactic radio signal 
expected from the strong intergalactic shocks in- 
volved in structure formation. We begin by con- 
structing the model from dimensional-analysis ar- 
guments [7J3] in H2.ll Next, we present a general- 
ized version of the model [HEj 1 corrected and cal- 
ibrated for a ACDM universe. The energy frac- 
tions £ e and £b are evaluated in H2.3I and a com- 
parison between model predictions and radio halo 
observations is presented in ij2.4l 

In the following we use a 'concordance' ACDM 
model: a flat universe with normalized vacuum 
energy density tt\ — 0.7, matter energy density 
VLm — 0.3, baryon energy density fls — 0.04, 
Hubble parameter h = 0.67, and an initial per- 
turbation spectrum of slope n = 1 and normal- 
ization erg = 0.9. The various parameters of the 
cosmological model are summarized in Tabled 
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Table 1 

Cosmological model parameters. 



Parameter 


Meaning 


Value 


h 


Hubble parameter 


0.67 


k 


Curvature 







Matter energy density 


0.30 




Dark matter energy density 


0.26 


n b 


Baryon energy density 


0.04 


n A 


Vacuum energy density 


0.70 


X 


Hydrogen mass fraction 


0.76 


n 


Fluctuation spectrum slope 


1 


OS 


Spectrum normalization 


0.9 



2.1. Dimensional Analysis 

Dimensional analysis arguments can be used to 
relate the mass M of a virialized halo, its velocity 
dispersion a, its smooth mass accretion rate M 
through strong shocks, its temperature T and the 
radius r s h of its strong accretion shock, by [TO] 



M = 



a/2 <t 3 (M,z) 
T~ GH{z) 



(1) 



M(M, z) = f a 



a 3 {M,z) 
3 G 



(2) 



2.5 x 10- 10 f acc h 70 a- 3 / 2 g(a)M yr" 1 , 



T(M,z) = frkg 1 f ia 2 {M,z) 



(3) 



1.8 x Wf T h% a 



2 /3 -1, 



{af' 3 Ml{ s K 



and 



r sh (M,z) = /, 



\/2 a(M,z) 
~5~ H{z) 



1.9f r h 7 Q /3 ag{a) 



-2/3 M l/3 



(4) 
Mpc . 



Here H(z) = 100 h a" 3 / 2 g{a) km s" 1 Mpc^ 1 is 
the Hubble parameter, with a = (1 + z) -1 and 

ff(o) = [n ro + n A a 3 + (i - n ro - o A )a] 1/2 , h 7Q = 

/i/0.70, /i ~ 0.65 m p is the average mass of a par- 
ticle (including electrons), M u = M/10 14 M Q , 
and k,B is the Boltzmann constant. The free pa- 
rameters face, /t and / r are dimensionless factors 
of order unity, which are roughly constant in the 
redshift range and cosmological model of interest, 



and must be calibrated separately. The defini- 
tions of the model parameters are summarized in 
Table along with their values as calibrated in 



Note that the above estimates refer to strong 
shocks only; weak shocks that may result from 
mergers of comparable mass objects are ignored 
since they lead to accelerated electron distribu- 
tions with little energy in highly relativistic elec- 
trons (for a discussion see [33]). Weak merger 
events may only enhance the predicted radio 
background from intergalactic shocks, and only 
in very low photon frequencies. 

Collisionless, non-relativistic shocks are known 
to accelerate electrons to highly relativistic en- 
ergies. In [35] it was shown, that SNR observa- 
tions suggest that the energy density of relativis- 
tic electrons accelerated by a structure formation 
shock constitutes a fraction £ e ~ 5% of the ther- 
mal energy density behind the shock (up to a fac- 
tor of ~ 2, see also discussion in W2. 3.1() . The 
maximal energy an electron can be accelerated 
to is limited by cooling, predominantly through 
invcrse-Compton scattering of background CMB 
photons, yielding a maximal electron Lorenz fac- 
tor lmax ~ 3 x 10 7 03| . 

The model focuses on strong shocks, which ac- 
celerate electrons to a power-law distribution of 
index p = 2 in the differential number of acceler- 
ated electrons per electron energy (equal energy 
per logarithmic interval of electron energy). The 
luminosity of a halo of mass M at redshift z due 
to inverse-Compton scattering of CMB photons 
may thus be estimated as 



uL?(M,z) 



2hi7„ 



n b M(M,z) 



£ e -k B T(M,z) 



1.1 x 10 42 I f acc f T 



0.05 



5/3 
l 70 



nit 



x a- 5 / 2 p(a) 5 / 3 l MU 3 



hi er g s 1 ■ ( 5 ) 



Assuming that the energy density of the down- 
stream magnetic field constitutes a fraction £b ~ 
1% of the downstream thermal energy density, the 
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magnetic field strength is given by 



B(M,z) = 0.14 



fr £,b 

/ 2 0.01 



1/2 



4/3 



'70 



7jh 

\l<rn 



1/2 



a 2 g 



(6) 



consistent with observation of galaxy cluster ha- 
los (see ^2.3.21 for discussion). This yields a syn- 
chrotron luminosity 

,L T (M, z) = B ^f^ ,Lf { M, z) 

1 9 X 10 39 ^ ^ acc ^T ^ e ^ B \ f, 13 / 3 



f? 0.05 0.01 



'70 



(7) 



7^ 



-5/2 



3(a) 



13/3 



,,7/3 _i 

MA erg s , 



where 



is the energy density of the CMB. 



Note that whereas the inverse- Compton signal 
depends on the parameter £g only logarithmi- 
cally (through the maximal energy attained by 
the relativistic electrons), the synchrotron emis- 
sion scales almost linearly with £b. 

Given the number density of halos of a given 
mass at a given redshift, dn(z)/dM, one may in- 
tegrate equation JSJ or equation J7J) to predict the 
inverse- Compton or the synchrotron background 
from intergalactic shocks, 



luLu) = / dz 



c dt 
dz 



dn(z) vL (M,z) 
dM 4tt(1 + z) 4 1 ' 



For example, with the Press-Schechter halo mass 
function and our cosmological model, we find 
an inverse-Compton background flux {vll 9) = 
1-3/occ/t (6/0.05) keV s" 1 cuT 2 sr" 1 , and a 
synchrotron background flux (vl^ yn ) = 3.2 x 
W~ 12 f acc f 2 fr 2 (6/0.05) (Cb/0.01) erg s" 1 cm" 2 ^- 

The model may be used to calculate the two- 
point correlation function of the radiation emit- 
ted by intergalactic shocks. The low optical depth 
of large clusters, which dominate the background, 
enables one to neglect cases where more than one 
cluster lies along the the line of sight, and approx- 
imate |70] 

J„(u)I„(v)) - (I v ) 2 

cdt f „ r dn(z) 1 



(9) 



dz 



dz 



dM- 



dM Trr sh (M,z)i 



vL v (M, z) 
47r(l + ^) 4 



Pi 



1|2 



ipd A {z) 



r sh (M, z) 



where u and v are unit vectors that satisfy 
u ■ v = cos"0. The function P^ 2 is the prob- 
ability that one line of sight passes through a 
halo of radios r s h , given that another line of sight 
passes through the same halo (see [3^] for detail). 
Halo-halo correlations, neglected in equation J?JJ, 
can only enhance the two-point correlation func- 
tion. The fractional correlation function, 6W = 
yj § 2 I V (ip) I {I v ) , is independent on the dynamical 
parameters /t, face, 6 and £s, and depends on 
the geometrical factor f r only through the rela- 

tion 2 U^fr) = fr'Ufr^A)- 

2.2. Fine- Tuned Model 

The Press-Schechter mass function is known to 
disagree with numerical simulations, predicting 
less rare, massive halos and more abundant, low 
mass halos. The Press-Schechter approach agrees 
better with cosmological simulations, if the naive 
critical over-density S c (z) indicating collapse (cal- 
culated for spherical collapse) is replaced by |5l)] 



5 c (z,M) = y/a~8 c (z)\l + b 



a 2 (M) 



a8 c (z) 



(10) 



where a(M) is the variance of the density field 
smoothed on a mass scale M, and a, b and c 
are parameters of order unity. Numerically, good 
agreement with ACDM simulations is obtained 
if one chooses a = 0.73, b = 0.34, and c = 0.81 
33 7 . Such a modification enhances the expected 
background from intergalactic shocks, because it 
increases the density of rare, massive halos, which 
dominate the extragalactic signal. We thus find a 
20% increase in the inverse Compton background, 
and a 50% enhancement in the synchrotron sig- 
nal. In what follows, we use the modified halo 
number density, unless otherwise stated. 

The spectrum of the integrated radiation from 
a halo accretion shock is essentially a broken 
power law, where the spectral break is introduced 
by inverse-Compton cooling. We focus on strong 
shocks which accelerate electrons to a power law 
energy distribution with index p = 2, implying 



2 At low [y < u br ) or high [y > u max ) frequencies (see 
'<\'2:2l . the scaling relation is more complicated. 
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that the photon spectrum is flat (vl v oc at 
high frequencies and scales as vl v oc v 1 ! 2 at low 
frequencies. The maximal Lorenz factor j max to 
which electrons are accelerated implies maximal 
synchrotron frequencies around 



,(M,z) 



3.6 x 10 1 



x(l 



B(M, z) 



-i 2 



0.1 fiG 
Hz . 



T(M,z) 
10 7 K 

(11) 



The spectral break frequency roughly corresponds 
to the minimal energy at which electrons manage 
to significantly cool (say, by rj = 50% of their 
initial energy) between the cosmic time of their 
host halo accretion shock (given by its redshift z) 
and the present epoch [55] . 



Vbr(M,z) 



9h, 



70 



n 



B(M, z) 
l-r/J 0.1/xG 

{1 + zf 2 

g( a ) 



(12) 



dz 



kHz 



The resulting synchrotron spectrum, integrated 
over all halos, cuts off below ~ 100 MHz. This 
implies that when including also weak shocks, 
the integrated spectrum vl v peaks at frequencies 
around 100 MHz. 

A hydrodynamical ACDM simulation |58I35| 
was used to obtain rough estimates of the dimen- 
sionless free parameters used in the model |36) . 
For this purpose, various global features of the 
simulated ACDM universe at the relevant epoch, 
< z < 2, were compared to their values accord- 
ing to the model. The calibration scheme was 
then tested by comparing features of the radia- 
tion from intergalactic shocks, as extracted from 
the simulation and as calculated from the model. 
The resulting parameter values are summarized 
in Table El The energy fractions £ e and £g can 
not be similarly evaluated from the simulation, 
but require independent observations (see H2.3( l. 

The temperature parameter fx was estimated 
using the temperature statistics of the baryonic 
component of the universe, indicating that al- 
though fx has a weak redshift dependence, it lies 
in the range 0.45 — 0.55 for the relevant epoch. 
The accretion rate parameter f acc was calibrated 
by the fraction of mass that has been processed 



by strong shocks in the relevant epoch. Compar- 
ing this fraction in the model and in the sim- 
ulation indicated that f acc has a weak redshift 
dependence, and lies in the range 0.08 — 0.17. 
The shock radius parameter f r has been esti- 
mated by studying the morphology of cluster ac- 
cretion shocks, suggesting that f r lies roughly in 
the range 0.6 — 1.2. 

An important aspect of the model is the sensi- 
tivity of the predicted synchrotron luminosity of 
a halo, but not its inverse-Compton luminosity, 
to halo asymmetry. Such asymmetries introduce 
spatial fluctuations in the thermal energy injec- 
tion rate through the halo accretion shock, and 
enhance the synchrotron emission which scales as 
the square of the thermal energy, vL^ 11 cx MT 2 . 
This effect can be incorporated into the model by 
introducing a geometrical correction factor |36| . 
The localized nature of hot regions along the 
shock, and the large ratio C — 10 between the 
temperature of these regions and the tempera- 
ture of the dimmer, more extended regions, jus- 
tify approximating all the emission from a halo 
accretion shock as originating from an effective 
smaller region, of scale r s h — ^s/i/C Replacing 
f r with a different parameter, f r = f r /Q, thus 
enhances the synchrotron luminosity of the halo 
by a factor £ 2 , while leaving its inverse-Compton 
luminosity unchanged. As an order of magnitude 
estimate, f r = f r /( ~ 0.05 - 0.20. Note that the 
two-point correlation functions are thus altered, 
power shifting to smaller angular scales. 



Table 2 

Intergalactic shock model parameters. 



Par. 


Quantity parameterized 


Value (Range) 


fr 


Temperature 


0.50 (0.45-0.55) 


face 


Accretion rate 


0.12 (0.08-0.17) 


fr 


Shock radius 


0.9 (0.6-1.2) 


fr 


Injection scale 


0.1 (0.05-0.2) 


& 


Electron energy fraction 


0.05 (0.02-0.10) 




Magnetic energy fraction 


0.01 (0.005-0.04) 



Dynamical parameters are defined in W2.l\ and 
calibrated in H2.2I Energy fractions are out of 
the shock thermal energy, and evaluated in 
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With the ACDM-modified halo mass func- 
tion and the best fit values described above, 
the model yields an inverse-Compton background 
vlf ~ 0.1(£ e /0.05) keV s" 1 cm" 2 sr" 1 . The 
synchrotron spectrum is presented in Figure 
the two-point correlation function in Figure ^ 
and the angular power spectrum (see H.'-t.l|l in Fig- 
ure 01 It is encouraging to note that after cali- 
brating the free parameters with essentially global 
features of the cosmological simulation, the model 
predictions agree well with results extracted from 
the simulation, regarding the 7-ray background 
|35| . the radio background, and the synchrotron 
two-point correlation function |37| . A more accu- 
rate calibration of the parameters, including an 
evaluation of their redshift dependence, may be 
obtained from a detailed analysis of the various 
clusters identified in cosmological simulations. 

2.3. Energy Conversion Efficiency 

Intergalactic shock waves are expected to accel- 
erate electrons to highly relativistic energies, and 
to strongly amplify magnetic fields. The average 
fractions of shock thermal energy deposited in rel- 
ativistic electrons (£ e ) and in magnetic fields (£b) 
are both important parameters of our model, each 
bearing linearly upon the predicted synchrotron 
signal, and thus deserve a special discussion. 

2.3.1. Electron acceleration 

Collisionless, non-relativistic shocks are known 
to Fermi accelerate power-law energy distribu- 
tions of relativistic particles. This phenomenon 
has been observed in astrophysical shock waves 
on various scales, such as in shocks forming when 
the supersonic solar wind collides with planetary 
magnetospheres, in shocks surrounding SNRs in 
the interstellar medium, and probably also in 
shocks in many of the most active extragalac- 
tic sources, quasars and radio galaxies 
The electron power-law distributions extend up 
to ~ 100 TeV energies in SNRs where shock 
velocities are of order v ~ 10 3 km s _1 , similar to 
intergalactic shock velocities. 

Although no existing model credibly calculates 
the acceleration efficiency, a simple argument |35| 
suggests that we may evaluate £ e in strong in- 
tergalactic shocks using the estimated accelera- 




te [deg] 



Figure 1. Fractional two-point correlation func- 
tion £u(if>, v) of synchrotron background from in- 
tergalactic shocks. The signal, calculated ac- 
cording to the calibrated model (see and 
Table 0) and shown for frequencies 1 MHz, 
100 MHz and 10 GHz (solid lines, bottom to 
top), scales roughly according to ^ v (ip, v ; f r ) = 

fr 1 £"(fr 1 i'> l ''i 1 ) for V br < V < V max , The 

signal calculated from a ACDM simulation (for 
v = 100 MHz, neglecting emission below the 
break frequency v\y ri circles with dashed line to 
guide the eye, error bars represent statistical er- 
rors introduced by averaging over random pairs 
of lines of sight |37]) is similar to the correspond- 
ing model result (dash-dotted line) , yielding even 
more small scale power. 



tion efficiency of SNR shocks. Consider an ide- 
ally strong, non-relativistic shock wave, such that 
the shock Mach number T>1 and the thermal 
energy of the upstream plasma is negligible with 
respect to the shock energy. The physics of such 
a shock is essentially determined by three param- 
eters: the shock velocity v, the upstream plasma 
number density n u , and the upstream magnetic 
field strength B u (in principle, the result may 
also depend on the detailed structure of the up- 
stream magnetic field). The upstream density 
may be eliminated from the problem altogether 
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by measuring time in units of v~ , where v p is 
the plasma frequency. The upstream magnetic 
field strength, parameterized by the cyclotron fre- 
quency, u c , can not be additionally scaled out of 
the problem. However, comparing v c .i (where 
subscript i denotes a property of the ions) to 
the growth rate of electromagnetic instabilities 
in the shocked plasma, Vi ns = v p ^v/c, indi- 
cates that their ratio in strong shocks satisfies 
{Vc,i/Vins? = {Bl/8ir)/(n u m p v 2 /2) < 1. We 
thus assume that there is a well behaved limit 
when this ratio approaches zero, implying that 
the upstream magnetic field has little effect on the 
behavior of strong shocks. With this assumption, 
we expect to find much similarity between suf- 
ficiently strong shocks in different environments, 
provided that their shock velocities are compa- 
rable, regardless of the plasma density and the 
strength or structure of the upstream magnetic 
field. The little effect of upstream magnetic fields 
on strong shocks is supported by recent SNR ob- 
servations, which suggest that the shocks produce 
magnetic fields much stronger than found far up- 
stream (see H2.3.2|l . 

We thus estimate the efficiency of strong in- 
tergalactic shocks from the strong SNR shocks, 
drawing upon the similarity between the veloc- 
ities of the two families of shocks. The energy 
fraction deposited in relativistic electrons by SNR 
shocks was estimated by several authors, the most 
reliable estimates found for remnants with 7-ray 
detection such as SNR1006 62 . From such stud- 
ies, we deduce an electron acceleration efficiency 
around £ e ~ 0.05, with an uncertainty factor of 
~ 2 (see jHSl for details). An independent, less re- 
liable method for estimating £ e relies on the ratio 
between the energies of cosmic-ray electrons and 
cosmic-ray ions in the interstellar medium, sug- 
gesting that £ e ~ 1% — 2% (for a discussion, see 
|35| ) . However, since the relation between this 
ratio in the interstellar medium and immediately 
behind intergalactic (or SNR) shocks is unknown, 
this estimate is highly uncertain. 

We stress that our estimate of £ e is based only 
on observations of SNRs, without provoking any 
elaborate model for the acceleration mechanism 
of the electrons by the shocks. The high efficiency 
of electron acceleration suggests that a non-linear 



theory (e.g. [HI and the references therein]) is re- 
quired in order to account for the shock structure 
and the particle acceleration process. Non-linear 
theories for diffusive acceleration of particles by 
shock waves are at the present stage incomplete 
(for a recent review, see |45| ) . In particular, such 
models fail to produce the flat electron spectra, 
which are implied by the observations and pre- 
dicted by a linear, test-particle approximation. 

2.3.2. Magnetic field amplification 

Magnetic fields in galaxy cluster halos have 
been estimated using several techniques, based 
on diffuse synchrotron emission from the cluster, 
preferably combined with inverse-Compton de- 
tection, on Faraday rotation of background or em- 
bedded polarized radio sources, and on the obser- 
vation of cold fronts in cluster X-ray images (for 
reviews, see 1^81311141 ) . The synchrotron emis- 
sion from a cluster measures the volume-averaged 
magnetic field weighted by the relativistic elec- 
tron distribution, and is thus most relevant for 
our model. Studies of the emission suggest vol- 
ume averaged magnetic field strengths of the or- 
der of a few 0.1 fiG, close to the values inferred 
for a minimal energy configuration with equal en- 
ergy in magnetic fields and in relativistic elec- 
trons. Faraday rotation and cold front studies 
suggest stronger magnetic fields of several /J.G, 
but are sensitive to different measures of the mag- 
netic field and to the assumptions made (see |36j 
for discussion). 

A fraction £g — 0.01 of a shock thermal en- 
ergy transferred to magnetic energy, implies mag- 
netic field strengths that are an order of magni- 
tude lower than their equipartition value (with 
respect to the thermal electrons, hereafter), and 
are close to their value in a configuration with 
equal energy in magnetic fields and in relativis- 
tic electrons (assuming £ e ~ 0.05). In clus- 
ter halos, our calibrated model reproduces, for 
this choice of £b, a volume-averaged magnetic 
field strength B ~ 0.1 ^G (for M ~ 10 14 M o ), 
although the 'bright spots' could contain mag- 
netic field strengths as high as 1 /iG. The low- 
est observational estimates for cluster magnetic 
field strengths are > 0.05 ^G, suggesting that £b 
could be higher than our chosen value by a factor 
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of a few, but is unlikely to be smaller than it by 
more than a factor of ~ 2. 

It is important to note that recent studies sug- 
gest that SNRs contain strong, near equipartition 
magnetic fields. For example, the highly local- 
ized nature of hard X-ray emission in the resolved 
sheets composing the shock front of SNR1006 
strongly suggests the presence of strong, near 
equipartition magnetic fields behind the shock 
(e.g. The narrow extent of emission ob- 

served upstream of these sheets implies a small 
upstream diffusion constant, corresponding to up- 
stream magnetic fields of strength > 10 /iG, with 
strong fluctuations (8B/B ~ 1) on small, d -C 
10 17 cm scales Such magnetic field fluctu- 
ations are far stronger, and have much smaller 
scales, than found in the surrounding ISM, and 
so must be induced by the shock. These conclu- 
sions strongly suggest that the assumption, that 
the physics of a strong shock depends only weakly 
upon the far upstream magnetic field, is valid, 
justifying the analogy between strong intergalac- 
tic shocks and strong SNR shocks of comparable 
velocities. 

2.4. Radio Halo Observations 

Radio halos are observed in ~ 35% of the 
young, massive galaxy clusters (with X-ray lumi- 
nosity L x > 10 45 erg s" 1 , see |24l22p . and their 
radio luminosity is known to be correlated with 
the cluster temperature 0U] an d mass - Fig- 
ure[21shows the specific luminosities of such halos 
for v = 1.4 GHz, plotted against the tempera- 
ture ranges of their host clusters. A simple power 
law fit to the data, L U (T, v = 1.4 GHz) cx 
gives 4> = 3.55, which is very similar to the power 
law index <f> — 7/2 predicted by our model, al- 
though the (calibrated) model yields a specific lu- 
minosity that is lower than observed by a factor 
of ~ 8 (equivalently, the temperatures are dis- 
crepant by a factor of ~ 1.8). A similarly good 
agreement between the model and the data is ob- 
tained when replacing the cluster temperatures 
with their estimated masses, giving L v (M 1 v = 
1.4 GHz) oc M 2 - 2 j2f2, similar to the model pre- 
diction L s J> n (M) cx M 7 / 3 . However, the esti- 
mated mass of a cluster is somewhat less certain 
than the measured temperature, and is sensitive 



to the definition of the cluster boundary. 
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Figure 2. Specific luminosities of radio halos for 
v = 1.4 GHz, plotted against the temperature 
ranges of their host galaxy clusters (error bars 
with cluster names, adopted from 26 221). The 
best power law fit to the data (solid line) scales 
as L„(T, v = 1.4 GHz) cx T 3 55 . The specific lu- 
minosity according to the model scales similarly, 
L V {T) cx T 7 / 2 , but the normalization (of the cali- 
brated model, dashed line) is different; the dash- 
dotted line presents the result of the model where 
the temperature was re-scaled by a factor of 1.8. 



Some comments regarding Figure are in 
place. First note, that the agreement between 
model and observations regarding the value of 
4> is insensitive to the free parameters of the 
model. The specific synchrotron luminosity of 
the model, when written as a function of halo 
temperature, scales linearly with the combina- 
tion faccfr 2 ! /t £e£s- Tne stron g dependence 
of this result upon the model parameters, and 
the low redshifts (z < 0.3) of all the halos shown 
in Figure [3 suggest that the discrepancy be- 
tween model and observations may be eliminated 
by better calibrating the parameters, and mod- 
elling their redshift dependence (note that this 
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will increase the radio signal we predict). Sec- 
ond, note that the low fraction of clusters with 
observed radio halos could result from a com- 
bination of two effects: (i) not all clusters have 
been significantly accreting new mass at the cos- 
mic time when they are observed; and (ii) the 
surface brightness of the halos observed is low, 
close to the instrumental sensitivity, suggesting 
that more radio halos will be observed with future 
telescopes. A varying accretion rate will tend to 
enhance the luminosity of the halos that are ob- 
served, possibly accounting for some of the dis- 
crepancy seen in Figure [3 Finally, note that in 
|46| . a cosmological simulation was used to show 
that (f> — 2.6 — 2.8 (for emission from electrons ac- 
celerated at shocks; for emission from secondary 
electrons produced by p-p collisions of cosmic ray 
ions, they found <fi — 4.1 — 4.2). However, the 
simulation of differs substantially from our 
model and simulation 35 37 , regarding their cos- 
mological model (SCDM), electron acceleration 
efficiency (£ e < 0.5%, see also 03), an d mag- 
netic field normalization ((B 2 ) 1 / 2 ~ 3 /iG for a 
Coma- like cluster). 

3. Observational Consequences 

In this section we examine the capability of 
present and next generation radio telescopes, in 
particular the SKA, to detect the radio emission 
from intergalactic shocks. In M3.ll we present the 
notations used. Various foreground and back- 
ground signals are discussed in ^3. 21 in an attempt 
to assess the optimal conditions for detecting the 
emission from intergalactic shocks. 

3.1. Notations 

As mentioned in S|21 the two-point correlation 
function of specific intensity fluctuations at a 
given angular separation ip, is defined as 

5 2 IM) = (6I v (n)SI v (v)) 

= (7„(u^(v))-^) 2 , (13) 

where u and v are unit vectors that satisfy u • v = 
cos ip, and SI v (u) = I v (u) — (/„). The specific in- 
tensity is related to the brightness temperature Xj, 



(which approximately equals the thermodynamic 
temperature for the sky brightness and the fre- 
quency range of interest) by l v = 2v 2 kBTb/c 2 . 
The following discussion of intensity fluctuations 
is thus equally applicable for temperature fluctu- 
ations, up to a multiplicative constant. 

It is often advantageous to study the Angular 
Power Spectrum (hereafter APS) of various sig- 
nals, in order to obtain a direct estimate of their 
importance at various angular scales. The an- 
gular power spectrum C; is defined through the 
relation 

oo 

5 2 IM) = T ^(2Z + l)GMP;(cos^) , (14) 
i=i 

where Pi (x) is the Legendre polynomial of degree 
I. The power at a given multipolc is often ex- 
pressed using its logarithmic contribution to the 
intensity variance, SIi = [1(21 + \)Ci/ Air] 1 / 2 (or 
STi for the variance of the brightness tempera- 
ture), where multipolc I roughly corresponds to 
angular scales 9 ~ 180°/Z. Using the orthogonal- 
ity of the Legendre polynomials, we may invert 
equation (JT3J| to find the APS, 

/■7T 

C 1 (v) = 2t o~ 2 I u (4))P l {cos4>)suvtl)dil) . (15) 
Jo 

The APS of synchrotron emission from intergalac- 
tic shocks may be calculated from the two-point 
correlation function (see Figure^) using equation 
i jTot . The resulting APS, presented in Figure 
peaks at multipoles / ~ 400 — 4000 for the rel- 
evant frequency range, corresponding to angles 
~ 3' - 30'. 

The sensitivity of a telescope is often estimated 
using its RMS noise level for a given beam size 
a = iriP 2 , denoted SI rms = [S 2 1(a)} 1 / 2 . For 
well behaved signals, where the two-point corre- 
lation function S 2 I(ip) is a monotonically decreas- 
ing function of ip, one finds 5 2 I(^p) < S 2 I rms (a). 
Hence, the RMS noise level imposes an upper 
limit to the two-point correlation function sen- 
sitivity, and may be compared directly with the 
logarithmic contribution to the variance at the 
corresponding multipole, SIi. 

3.2. Competing Signals 
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Figure 3. Angular power spectrum of the syn- 
chrotron background from intergalactic shocks, 
calculated according to the calibrated model 
(see and Table |2J for frequencies 1 MHz, 

100 MHz, and 10 GHz (dashed lines, bottom to 
top). The result of a ACDM simulation for v = 
100 MHz (neglecting emission below the spectral 
break |37], shown as a band with thickness cor- 
responding to the statistical error) exhibits more 
power than the corresponding result of the model 
(dash dotted line), in particular on small, ar- 
cminute scales. 



In the following, we examine various fore- 
ground and background signals in the low fre- 
quency radio band. We focus on Galactic 
synchrotron foreground and on discrete radio 
sources, which contaminate the expected fluctu- 
ating signal on large and on small angular scales, 
respectively. Other extragalactic low frequency 
signals — bremsstrahlung from Lya clouds and 
21 cm tomography — are also discussed. The 
contribution of each signal to the sky brightness 
is shown in Figure 01 and its contribution to the 
variance on 0°.5 angular scales in Figure 

3.2.1. Galactic Synchrotron Emission 



At low frequencies, v < 1 GHz, the brightness 
of the sky is dominated by Galactic synchrotron 
emission, produced by cosmic-ray electrons gy- 
rating in the magnetic fields of the interstellar 
medium (ISM). At low Galactic latitudes, there 
is also a substantial contribution of free- free emis- 
sion from low-latitude H II regions (e.g .1 ). In 
the frequency range of interest, v <§; 10 GHz, 
the two-point correlation function is dominated 
by Galactic synchrotron emission, because of its 
significant power on large angular scales (see [3U 
for discussion). We thus focus on the contami- 
nation by Galactic synchrotron emission in high 
Galactic latitudes and on small angular scales. 

With large uncertainties regarding the spatial 
distribution of Galactic cosmic rays and mag- 
netic fields, the best estimates of Galactic syn- 
chrotron emission are obtained by extrapolating 
direct measurements carried out at frequencies 
and angular scales where synchrotron emission 
is reliably measured. We thus make the key 
assumption, often used in the literature, that 
the multipolc dependence of the Galactic syn- 
chrotron APS varies little with frequency, such 
that Ci(v) ~ f{l)g(v). The frequency depen- 
dence g{v) of synchrotron emission can be mod- 
elled as a broken power law, and is discussed in 
|3f>| . The APS multipole dependence f(l) may be 
extracted from high resolution maps at ~ GHz 
frequencies down to arc-minute scales (multipoles 
I ~ 6000, [OBI)) an d ma y m principle be extrapo- 
lated into the entire frequency regime and multi- 
pole range of interest. 

The angular power spectra extracted from ra- 
dio maps depend on a number of factors, most 
importantly the Galactic latitude range examined 
and the efficiency at which discrete radio sources 
are removed from the map. Unfortunately, high 
resolution surveys are available mostly for low 
latitudes, which are less promising as potential 
search sites for an extragalactic signal because 
they exhibit a stronger Galactic foreground, in 
particular at small angular scales. Nonetheless, 
improved high and medium latitude data have re- 
cently made it possible to extend our knowledge 
of the high latitude APS up to multipoles I ~ 800. 

The high latitude (60.5° < b < 84.5°) low reso- 
lution (FWHM 0.6° - 2.3°) maps of 11 for 5 fre- 
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quencies between 408 MHz and 1411 MHz, have 
been used to estimate the APS up to I ~ 70 (e.g. 
[T31 table 6]). Analysis of the Bonn 408 MHz full- 
sky survey j3U] suggests a power-law APS of the 
form Ci(l > 1) ex l~P with ~ 2.5 - 3.0 down to 
the survey resolution limit 0.85°, corresponding 
to multipoles I ~ 200 63 . With a characteristic 
coherence scale ~ 5° — 10° for synchrotron emis- 
sion (e.g. |59I6| 1. Ci cx (I + 5)~ 3 was claimed 
to describe the APS rather well. Thus, al- 
though the contribution of Galactic synchrotron 
emission to the two-point correlation function is 
significant, it is dominated by large angular scales 
~ 5° - 10° (low multipoles I < 40) and intro- 
duces little contamination at small angular scales 
(see Figure 01 . The APS exhibits strong fluctu- 
ations across the sky, suggesting that 'quiet' re- 
gions may be identified and selected as preferable 
search sites for an extragalactic signal. 

An analysis of the 2.3 GHz Rhodes map |34| . 
with FWHM resolution 20', reveals strong vari- 
ations of the APS with Galactic latitude, and 
a large contribution from discrete radio sources 
which tend to flatten the APS At high lati- 
tudes (|6| > 20°), removal of point sources using 
median filtering gives 123 

Q ~ A 2 l~P , (16) 

with A = 0.3±0.2 K, and (3 = 2.92±0.07, valid up 
to I ~ 100. Although the quoted uncertainties are 
large, comparison with the low resolution studies 
mentioned above indicates that the normalization 
A is probably no larger than 0.3 K. Similarly 
steep APS, with (3 varying in the range 2.60 — 3.35 
for different medium latitude (|&| < 20°) regions, 
were found from a 1.4 GHz survey carried out 
with the Effclsberg 100 m telescope j^H], after re- 
moving > 5<7 unresolved point sources pQ. The 
high angular resolution of the Effelsberg survey, 
~ 9.35', thus suggests that the steep power law 
discussed above holds up to multipoles I ~ 800. 
Note that after removing the brightest intensity 
peaks associated with low latitude H II regions, 
the APS highly resembles the APS of the polar- 
ized intensity component pQ, as expected for syn- 
chrotron emission. 



3.2.2. Discrete Radio Sources 

Discrete radio sources, mostly radio galaxies, 
active galactic nuclei (AGNs) and normal galax- 
ies, make an important contribution to the ex- 
tragalactic radio background. Catalogues of ra- 
dio sources at several frequencies have been used 
to estimate their contribution to the radio sky 
[57173] . This procedure is limited by the un- 
certain contribution of faint, unidentified sources 
(mostly normal galaxies), which dominate the 
background (see for example [73]) and must be 
modelled. Other studies have estimated the con- 
tribution of discrete sources to the radio back- 
ground, utilizing the well-known correlation be- 
tween the infra-red and radio flux densities of in- 
dividual galaxies 50 28 . However, the results of 
such studies are sensitive to the unknown redshift 
evolution of the sources. In addition, the radio- 
infra-red correlation holds only in cases where the 
radio emission is associated with star formation. 
Hence, such studies must be supplemented by 
independent estimates of sources not associated 
with star formation, such as a catalogue-based 
estimate of the emission from AGNs |53| . 

The studies outlined above suggest that the 
contribution of discrete sources to the radio sky 
is roughly an order of magnitude lower than 
the observed sky brightness. This corresponds 
to 25% — 65% of the total extragalactic back- 
ground, as estimated by spectral modelling of low 
frequency radio observations |15I13I36| . although 
some evolutionary models (e.g. (SU|) may even 
account for the entire strong extragalactic signal 
calculated by ^3 • The spectral index of the inte- 
grated emission from point sources has been esti- 
mated to lie in the range s = 0.7 — 0.8 for frequen- 
cies in the range 100 MHz - few GHz |57I39I28| . 
although some dependence upon scale may be ex- 
pected because the source distribution is essen- 
tially bimodal [Sj. At very low [y < 10 MHz) 
frequencies, the background from discrete sources 
turns around, probably because of synchrotron 
self-absorption |57j . 

It is easier to estimate the contribution of point 
sources to the APS, because the latter is domi- 
nated by the bright, well studied sources. Dis- 
crete radio sources have approximately a Pois- 
son distribution in the sky, because projection 
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Figure 4. The low frequency radio sky. Ground-based (circles and error bars for the Galactic po- 
lar regions) and space-based (triangles for IMP-6 spacecraft minimum and maximum) observations are 
dominated by Galactic synchrotron emission (dash dotted line shows the modelled Galactic foreground 
toward the polar regions jSE]) m low frequencies, and by the CMB (double-dashed line) at high fre- 
quencies. Various estimates of the extragalactic component (legend, see text) are typically an order of 
magnitude lower than the Galactic foreground. The background from intergalactic shocks according to 
the calibrated model (of H2.2l and Tabled solid line) is roughly of the same magnitude, scaling according 
to I, y cx faccfxf^ 2 - This signal is also shown according to a ACDM simulation (accounting only for 
emission above the spectral break, dotted line, [HZ])- The frequency ranges of the SKA and the LOFAR 
(double arrows and corresponding shaded regions) and brightness temperature contours (dotted lines, 
labelled on the top and on the right axes) are also shown. 
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through their wide redshift distribution effec- 
tively diminishes their correlations (e.g. 63 ). 
Hence, as long as the angular scales concerned 
are much larger than the angular extent of the 
sources, the two-point correlation function van- 
ishes and the APS is flat, C; oc I (white noise), 
such that roughly SIi oc I. One may estimate the 
APS as (e.g. 

ci= L°° s2d 4 ds > (i7) 

where S is the source flux and dN / OS is the dif- 
ferential number density of point sources in the 
sky. The limited sensitivity S m i n of the observa- 
tions limits our knowledge of dN/dS to sources 
with flux S > S m in, but since the number den- 
sity of faint sources is not too steep, the possi- 
ble error introduced is small. The upper limit of 
the integral, essentially determined by the bright- 
est sources, can be lowered in order to reduce 
the noise, by modelling and removing the bright- 
est sources from the analyzed map. However, 
this quickly becomes laborious as the number of 
sources increases, and introduces inevitable errors 
associated with source removal uncertainties. 

Following [fi3| . we use source counts produced 
by the 1.4 GHz VLA FIRST all-sky survey. A 
small fraction (~ 10~ 4 ) of the discrete sources in 
the FIRST catalogue have angular scales larger 
than an arc-minute |71| . so the use of equation 
1|17JI is legitimate when dealing with multipoles 
' IS 10 4 - The number density of radio sources 
(after performing a resolution correction for ex- 
tended sources, [7T]) is well fit in the flux range 
5=1 mJy — 1 Jy by 

For faint sources in the range 5=10 /iJy— 1 mJy, 
the number count steepens to dN/dS oc 5" 7 with 
7 = 2.2 ± 0.2 74 . This changes the integrated 
brightness considerably, but has a negligible effect 
on the APS. Using equation l|17|l. we find that for 
a modest cut at S max = 100 mJy, correspond- 
ing to removal of the ~ 7 x 10 4 brightest sources 



all-sky, Ci{v) ~ 1.7 x 10~ & (v/lA GHz)- 2 75 K 2 . 
Using the 6C survey carried out at 151 MHz |2~9"] . 
we find that the above cut is equivalent to removal 
of sources brighter than ~ 3 Jy at 150 MHz. As 
shown in Figure |SJ such a cut places the emis- 
sion from point sources at the same level as the 
Galactic foreground for angular scales ~ 0°.5. At 
smaller angular scales, contamination from dis- 
crete sources becomes increasingly worse, requir- 
ing the removal of more sources. 

3.2.3. Bremsstrahlung from Lya Clouds 

The thermal bremsstrahlung emission from 
Lya clouds has been calculated in These 
clouds are most likely in photo-ionization equi- 
librium with the UV background, permitting an 
estimate of the bremsstrahlung signal if the UV 
background intensity is known. UV background 
intensity estimates based on the proximity effect 
|55| imply a bremsstrahlung signal |36| 

J v = 10" 22±1 erg s" 1 cm" 2 sr" 1 Hz" 1 , (19) 

for v <C 10 14 Hz. This result neglects the contri- 
bution of optically thick clouds at high (z > 5) 
redshifts, and possible additional ionization fields, 
such that the actual bremsstrahlung flux may be 
somewhat higher than given above. 

The amplitude of intensity fluctuations intro- 
duced by the finite number of H I clouds along 
any given line of sight has also been evaluated 
|41| . Assuming a random cloud distribution in 
redshift and in column density, recent UV back- 
ground estimates imply RMS fluctuations 36 

AJ„ = 10" 23±1 erg s" 1 cm" 2 sr" 1 Hz" 1 , (20) 

for v <C 10 14 Hz. Such fluctuations are expected 
on ~ 7" scales (corresponding to I ~ 10 5 ), al- 
though weaker fluctuations introduced by clouds 
associated with large-scale structure may appear 
even on ~ 10' scales. This fluctuation level should 
be regarded as a lower limit, because in addition 
to the assumptions leading to equation l|19|l . the 
non-randomness of the cloud distribution further 
enhances the fluctuation signal. Nevertheless, 
without a unique spectral or temporal signature, 
confusion with the other signals discussed in this 
section, in particular discrete sources, is likely to 
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Figure 5. Logarithmic contribution SIi of various signals to the variance at multipolc I = 400 (cor- 
responding to 6 ~ 0°.5). The background from intergalactic shocks is shown according to the cali- 
brated model of H2.2l and Tabled (solid line), and scales roughly according to SIi(l; £ e , £b, face, fr, fr) = 
^e^Bfaccfrfr 3 &h{frh 1, 1, 1, 1, 1)- The signal according to a ACDM simulation (diamonds with statis- 
tical error bars for I = 400, 10 3 and 10 4 , from bottom to top, accounting only for emission above the 
spectral break, see |HZJ) scales linearly with the combination £ e £b- The high- latitude Galactic foreground 
(dash-dotted line) and the integrated emission from discrete radio sources (dotted line, assuming removal 
of sources above S cu t = 100 mJy at 1.4 GHz) are presented. Above ~ 10 GHz, CMB fluctuations take 
over (dashed line extrapolated from WMAP measurements, error bars present detections by the Cam- 
bridge Radio Telescope on 0°.2 — 0°.5 scales, and by the Owens Valley Radio Observatory on 0°.l — 0°.6 
scales; for references see |H2I72| V Bremsstrahlung emission from Lya clouds (low horizontal shaded re- 
gion) and the maximal 21 cm tomography signal (dashed contour) are seen to be much weaker than 
the above competing signals. Continuum surface brightness sensitivities of the SKA (square) and the 
LOFAR (star), and the side-lobe confusion limit of the VLA (triangles, for 10 minutes in configuration 
D, see http://www.vla.nrao.edu/astro/guides/vlas) are shown (labels denote beam widths). Brightness 
temperature contours are shown as dotted lines (labelled on the top and the right axes). Emission from 
intergalactic shocks is seen to dominate the sky at this angular scale for frequencies v < 500 MHz, pro- 
vided that £ e £b ^ 10~ 4 . According to the simulation, the emission dominates the sky at this frequency 
range for 9 < 0°.5, even if £ e £ B = 5 x 10~ 5 . 
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preclude detection of the bremsstrahlung signal 
in the relevant frequency range. 

3.2.4. IGM 21 cm Tomography 

The 21 cm (1.4 GHz) spin flip transition of 
atomic hydrogen was proposed as a probe of the 
IGM at the epoch prior to reionization jll] (see 
also jHTj) . Spatial inhomogeneities in the IGM 
may be observed today as rcdshiftcd emission or 
absorption fluctuations against the CMB, in situ- 
ations where early sources of radiation decoupled 
the IGM spin temperature from the CMB tem- 
perature. In principle, the combined angular and 
spectral signal, stronger than the CMB fluctua- 
tions by two orders of magnitude, may be used 
to trace the 'cosmic- web' structure of the early 
universe, in both space and cosmic time. 

However, as pointed out by JS|> at the rel- 
evant frequency range 50 — 200 MHz, contam- 
ination by radio point sources imposes a seri- 
ous contamination even for the maximal signal 
amplitude AT ~ 10 mK expected. A spectral 
~ 10 mK feature, caused by the fast rise of the 
Lya background as the first UV sources reion- 
ized the IGM, could possibly be detected behind 
the spectrally continuous foreground |75I48I42_2*51 
and references therein]. Emission from inter- 
galactic shocks provides an additional important 
source of confusion for 21 cm tomography. 

4. Discussion 

We have calculated the imprint of intergalac- 
tic shocks on the low frequency radio sky, us- 
ing a simple model, calibrated and tested with 
a ACDM simulation. Our main results are il- 
lustrated in Figures 0] and [5] (for the calibrated 
model parameters summarized in Table EJ). Fig- 
ure 01 shows the contribution of various signals 
to the low frequency radio sky, suggesting that 
emission from intergalactic shocks contributes up 
to a few tens of percent of the extragalactic radio 
background below 500 MHz. Figure [5] depicts the 
angular power spectrum of various signals on an 
angular scale of ~ 0°.5, along with the sensitivity 
and the angular resolution of the SKA. 

We conclude that the design of the SKA, as 



well as the design of next generation radio tele- 
scopes such as the LOw Frequency Array (LO- 
FAR 3 ) and the Astronomical Low Frequency Ar- 
ray (ALFA 4 ), are more than sufficient for detec- 
tion of the angular fluctuations introduced by in- 
tergalactic shocks, as calculated in (see also 
36 37 ). Identification of the signal is limited by 
confusion with Galactic foreground and with dis- 
crete radio sources. Foreground fluctuations in 
the synchrotron emission of our Galaxy constrain 
a clear detection of the signal to sub-degree scales 
and to high Galactic latitudes. Confusion with 
discrete radio sources requires that the brightest 
sources be modelled and removed. With a feasible 
point source cut (100 mJy at 1.4 GHz, or 3 Jy at 
150 MHz), the signal dominates over the compet- 
ing signals at angular scales 10' < 9 < 1° (in low 
frequencies, for £ e £.B > 5 x 10~ 5 , see discussion 
below), whereas detection on arcminute scales 
will require a more ambitious point source re- 
moval. The spectrum of intergalactic shock emis- 
sion indicates that the signal is most pronounced 
at frequencies around ~ 100 MHz, planned to be 
covered by the SKA, the LOFAR and the ALFA. 

The calculated signal is sensitive to uncertain- 
ties in the model parameter calibration. The log- 
arithmic contribution to the variance (shown in 
Figure EJ) scales roughly according to 

_e, £b i face, h, Jr) (21) 

= Us^^Shifrl; 1,1,1,1,1). 
fr 

The radio signal calculated from the cosmolog- 
ical simulation |37], however, depends only on 
the energy fractions £ e and £g, scaling linearly 
with the combination £ e £g . Comparison between 
the predictions of the model and of the simula- 
tion (e.g. in Figures ^ El - EJ suggests that C; 
was in fact underestimated by our parameter cal- 
ibration scheme. This implies, for example, that 
intergalactic shocks introduce intensity fluctua- 
tions of magnitude SIi > 8 x 10 _19 (^e€s/5 X 
10- 4 )(i7l00 MHz)" 1 erg s" 1 cm" 2 sr" 1 Hz _1 on 
multipoles 400 < I < 2000, corresponding to 
temperature fluctuations STi > 260(£ e £s/5 x 

3 see http:/ /www.lofar .org| 

4 see http:/ /sgra.jpl. nasa.gov/html_dj/ALFA. html 
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10" 4 )(^/100 MHz)" 3 K. 

The fluctuations introduced by intergalactic 
shocks, as estimated above, dominate the sky at 

< 500 MHz frequencies if £, e £,B ^ 10 -4 , and are 
dominant even at frequencies as high as 10 GHz 
if £ e £s 3 x 10~ 4 . The simulation suggests 
that the signal is stronger than predicted by 
the model, in particular on small angular scales. 
Therefore, according to the simulation, emission 
from intergalactic shocks will dominate the sky at 

< 0°.5 scales and < 500 MHz frequencies, even if 
^eCs — 5 x 10~ 5 . We have used observations of 
SNR shocks and of magnetic fields in the halos of 
galaxy clusters, in order to show that for strong 
intergalactic shocks £ e £g — 5 x 10~ 4 , unlikely to 
be smaller than this value by more than a fac- 
tor of ~ 4 (see i!2.3|) . Detection of the signal is 
not unrealistic even in case 5Ii has been mildly 
overestimated. A signal lower than calculated in 
S0 by a factor of a few, may still be identified 
at ~ 1' — 10' scales, if faint point sources are 
modelled and removed. If the signal was overes- 
timated by ~ 2 orders of magnitude, it may still 
be detectable by next generation radio telescopes 
such as the SKA, by means of cross-correlation 
with known tracers of large-scale structure, such 
as galaxy counts, and in the future, with 7-ray 
emission from intergalactic shocks. 

Interestingly, detection of the intergalactic 
shock signal is just possible with present high res- 
olution radio telescopes, such as the VLA (at the 
maximal sensitivity configuration, see Figure [5j. 
The calculated signal should be detectable at high 
latitudes, ~ 100 MHz frequencies and sub-degree 
scales. The signal may have already been de- 
tected by CMB anisotropy studies, at ~ 10 GHz 
frequencies and 9 < 1° scales [70], and in quiet re- 
gions of ~ GHz high resolution surveys at 1' — 10' 
scales [30]. It will probably be easiest to iden- 
tify the signal at arcminute scales, by modelling 
point sources and cross-correlating the maps with 
known tracers of large-scale structure. 

Future detection of radio emission from inter- 
galactic shocks will have important implications 
on our understanding of cosmology and astro- 
physics. Detection of the signal will provide the 
first identification of intergalactic shocks, reveal- 
ing the underlying cosmological flows and provid- 



ing a test for structure formation models. The 
signal, in particular when combined with 7-ray 
detection, will provide a measure of the unknown 
magnetic fields in the intergalactic medium. Al- 
though non-trivial for interpretation, such a mea- 
sure of the magnetic field may provide insight into 
the unknown processes leading to IGM magneti- 
zation. The signal may confirm the existence of 
the undetected warm-hot intergalactic medium, 
and provide a tracer for its distribution. 

Synchrotron emission from intergalactic shocks 
is correlated with the large-scale structure of 
the low-redshift (z < 1) universe, tracing young 
galaxy clusters and filaments. The signal could 
thus account for some observed features of the 
radio signature of galaxy clusters, namely ra- 
dio halos and radio relics. Extended accretion 
shocks could contribute to radio halos (e.g. to 
the large radio halo of the Coma super-cluster, 
see and the references therein]), whereas lo- 
calized shocks (e.g. where galaxy filaments chan- 
nel large amounts of gas into the cluster regions) 
may be responsible for some radio relics observed 
at the outskirts of clusters (such as the prototype 
relic found in the Coma super-cluster, 1253+275 
|20p. and perhaps also for some of the anoma- 
lous features observed in radio halos (e.g. in the 
unrclaxed clusters described in |27p. 

Radio emission from intergalactic shocks is an 
important source of contamination for other radio 
signals, such as the low frequency CMB, Galac- 
tic synchrotron fluctuations on sub-degree scales, 
and competing extragalactic radio signals such 
as bremsstrahlung from Lya clouds, and 21 cm 
tomography. Emission from intergalactic shocks 
should be taken into account when evaluating the 
propagation of ultra-high energy photons with en- 
ergies above 10 19 eV, where the effect of the radio 
background on the transparency of the universe 
is stronger than the effect of the CMB. 
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